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Abstract We discuss the regularity of solutions of 2D incompressible Navier- 
Stokes equations forced by singular forces. The problem is motivated by the 
study of complex fluids modeled by the Navier-Stokes equations coupled to a 
nonlinear Fokker-Planck equation describing microscopic corpora embedded 
in the fluid. This leads naturally to bounded added stress and hence to 
W~ 1,oc forcing of the Navier-Stokes equations. 
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1 Introduction 

We discuss the regularity of solutions of 2D incompressible Navier-Stokes 
equations forced by singular forces. The problem is motivated by the study 
of complex fluids modeled by the Navier-Stokes equations coupled to a non- 
linear Fokker-Planck equation describing microscopic corpora embedded in 
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the fluid. This leads naturally to bounded added stress and hence to W~ 1,co 
forcing of the Navier-Stokes equations. A more detailed description of the 
problem in question, together with an application of the results in the present 
paper can be found in our forthcoming paper [2]. 

In this paper we focus on the 2D Navier-Stokes issues. The global ex- 
istence of energy solutions and their uniqueness are well known as classical 
results of J. Leray for the Cauchy problem and O. Ladyzhenskaya for initial 
boundary value problems in bounded domains. These results remain to be 
true for singular forces as well. 

The regularity of energy solutions with relatively smooth forces is also 
known. Regularity can be established, for instance, by scalar multiplication 
of the Navier-Stokes equation by the Stokes operator of the velocity field, 
integration by parts, and application of Ladyzhenskay's inequality 

IMIi*(R2) < v / 2||m||l2 (M 2 ) ||Vm|| L 2 (R 2 ) , Vm G C^(IR 2 ). 

This procedure yields summability of the second spatial derivatives. Further 
regularity can be obtained perturbatively, with the help of the linear theory. 

The regularity of energy solutions with singular forcing is limited. The 
best one can expect is Holder continuity of the velocity field. We prove Holder 
continuity at a local level, in both space and in time. We assume that our 
local solution has finite energy and the pressure field is in L 2 . This latter 
assumption seems restrictive: we are not able to justify it for general initial 
boundary value problems with reasonable singular forcing. The assumption 
is however satisfied in the absence of boundaries, i.e., for the Cauchy problem 
in the whole space and for the initial value problem on the torus. We briefly 
explain in this paper how the local regularity results can be applied to the 
Cauchy problem in the whole space. 

In our proof, the Holder continuity of the velocity field depends quanti- 
tatively on the modulus of continuity of the function uj i— ► J \u\ 4 dz. In order 

to be able to apply this regularity result to coupled systems or to families of 
Navier-Stokes systems, this modulus of continuity needs to be a priori uni- 
formly controlled. We achieve this in the absence of boundaries by obtaining 
higher integrability of the velocity, u G L°° (dt; L r (dx)) , r > 4. In order 
to obtain the higher integrability we prove the generalized Ladyzhenskaya 
inequality that reads 

T 

|M|x2r (K 2) < — ^|k||^(]R2)||Vw||i2( R 2), Vw G C£°(IR 2 ) 
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for r > 2. The proof is elementary and can be found in the Appendix. 



2 Notation and Local Regularity Result 

We assume that Q and Qi are domains in IR 2 such that Q 1 <s Q and < 7\ < 
T, and let 

Q = Q x (-T, 0), Qi = fti x (-7\, 0). 

Parabolic balls will be denoted as Q(z , R) = B(x , R) x (t — R 2 , t ), where 
Zo — (%o,to), G I^ 2 , £ ^, and B(x 0: R) is an open disk in M 2 having 
radius R and centered at the point x . 
We use the following notation for mean values: 

Q(z ,R) B(x ,R) 

L P (Q) and W l,p (Q) stand for usual Lebesgues and Sobolev spaces of functions 
defined f2, and the norm of the Lebesgues space is denoted by || ■ || mi n. For the 
forcing we are going use a functional space M 2i7 (Q) with parameter < 7 < 1 
and seminorm 

1 f ~ 

M2MQ) = SU P Rl ~A \n(, mi / \f(. z )-U)z a ,R\ 2 dz) 2 < 00. 
Q{z ,R)cQ Mv(2 ,-Kj| J ' 

Q(zo,R) 

We denote by c all positive universal constants. Our regularity result can be 
formulated as follows. 

Theorem 2.1. Assume that we are given functions 

ueL 4 (Q;R 2 ), P EL 2 (Q), F G M 2j7 (Q; M 2x2 ) (2.1) 

with < 7 < 1, satisfying the Navier- Stokes equations 

d t u + u- Vu- Au + Vp = -divF, divw = (2.2) 

in Q in the sense of distributions. 
Then 

u E CPQQJ (2.3) 

if < 7 < 1 and 

u e BMO(Qi) (2.4) 

if 1 = 0. 
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Remark 2.2. The Holder continuity and the BMO space are defined with 
respect to parabolic metrics. 

Remark 2.3. The corresponding norms are estimated in terms of the quan- 
tities |M| 4 ,q, \\ph,Q, \\F\\m 2 ^(Q), dist(fii,<9fi), T-Ti, and the modulus of 
continuity of the function to i— > J \u\ A dz. 

Several additional results can be proved by means of minor modifications 
of the proof of Theorem 12. II Before stating one them, we define usual energy 
spaces for the 2D Navier-Stokes equations. Let H and V be completions of 
the set of all divergence-free vector fields from C^°(1R 2 ; 1R 2 ) with respect to 
the L 2 norm and the Dirichlet integral, respectively. 

Proposition 2.4. Let u G L°°(0, T; H) H L 2 (0, T; V), p E L 2 (0, T; L 2 (R 2 )) 
be a solution of the Cauchy problem 

d t u + u ■ Vu — Au + Vp = — div F, div u — 0, (2.5) 

«(.,0) = o(-)Gff, (2.6) 

where F 6 L q (Q T ; M 2x2 ) n L 2 (Q T ; M 2x2 ) with q > 4 andQ T = R 2 x (0, T) . 

Then, given < s < T , there exists a constant C depending only on s, 
the norms of F in L q (QT] M 2x2 ) and in L 2 (Qt]M. 2x2 ) and the norm of a in 
H such that 

|MU°°(R 2 x(s,T)) < C. (2.7) 

Moreover, the function u is Holder continuous in M 2 x [s, T] with exponent 
7 = 1-4. 

Remark 2.5. The existence and uniqueness of a solution to the Cauchy 
problems A2.5\) and A2.6\) with above properties is well known, see f^j. 

Remark 2.6. The same statement is valid in the case of periodic boundary 
conditions. More generally, it is true as long as the pressure field is in L 2 . 



3 Proof of Theorem 2.1 



We are going to analyze differentiability properties of the velocity field u in 
terms of the following functionals: 



$(«; zq, g) 



u - (u) Z0 J 4 dz) , ^(u;z ,q) 



\u\ dz 



Q(zo,q) 



Q(zq,q) 
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D(p;z ,g)= J \p- \p) Xo J 2 dz. 



Q(zo,e) 

The following two statements are well-known. 
Lemma 3.1. Let the function v G L 4 (Q(zq, R)) satisfy the heat equation 

d t v — Av = 

in Q(zo, R). Then 

^(v;z ,g)<c^y^(v;z ,R) (3.1) 

for all < g < R. 

Lemma 3.2. Given G G L 2 (Q(zq, R); M 2x2 ), there exists a unique function 
w G C([t - R 2 , t ]-L 2 {B{x Ql R)- IR 2 )) n L 2 ([t - R 2 , t ]; W l > 2 {B{x Ql R); K 2 )) 
such that 

d t w — Aw = — div G 

in Q(z , R) and 

w = 

on the parabolic boundary oJQ(zq, R). Moreover, the function w satisfies the 
estimates: 

1^11,0(^0,^) = SU P IK'^llLwi) + W^ W \\Iq(z ,R) 
t -R 2 <t<t 

< 2 l|G'll2,Q( Z o,-R)' ( 3 ' 2 ) 

<$>{w;z Q ,R) < c\w\l Q{z0tR) . (3.3) 

The next couple of statements are about some properties of the solutions 
of the system 

Lemma 3.3. Under the assumptions of Theorem \2.1\ we have 
$(u; zq, g) < c{ ^ 4 + ^{u;z ,R) <5>{u; z , R) + 

+ D(p;z ,R) + MR 2+2 ^ (3.4) 
whenever Q(zo, R) C Q and < g < R. Here, M = \\F\\m 2 tQy 
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Proof. Setting 

G = F - {F) Z0>R + (p - \pU, R )I + (u- (u) Z0 , r ) <g> it 
in Lemma 13.21 we get the following estimate for w 

<c J [\F- (F) Z0tR \ 2 + \p- \p] X0 , R \ 2 + \u- (u) Z(hR \ 2 \u\ 2 



dz < 



Q(z ,R) 



< c 



MR 2+2 ^ + D{p; Zq, R) + V(u; z , R)$(u; z , R) 



(3.5) 



Obviously, the function v = u — w satisfies the heat equation. Then, applying 
Lemma [3. 1[ we find 

$(u - w; z , g) < c(-|) $(m - w; z , R). 



The latter inequality gives us: 



zq, q) < c ( — ) $(u; z , + zb, #) 



(3.6) 



Combining (13.51) and f)3.6p . we arrive at f)3.4p and thus Lemma is proved. 
Lemma 3.4. Under the assumptions of Theorem \2.1\ we have the estimate 

D(p;z Q ,Q) < c[^yD(p;z ,R) + y(u;zo,RMu;z ,R) + MR 2+2 ^ (3.7) 

whenever Q(zq, R) C Q and < g < R. 

Proof. The crucial part of the proof is the pressure decomposition 

V = P1+P2, 

where the first component pi satisfies the identity 

J pi&ydx = - J - {u) Zo>R ) ®u + F - {F) ZQ> R^ : V 2 (pdx 

B{x ,R) B(x ,R) 
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where the test function ip e W 2 2 (B(x , R)) is subject to the Dirichlet bound- 
ary condition: ip = on 8B[xq, R). It is not difficult to show that such a 
function p\ exists and obeys the estimate 



J \Pl-\Pl]xo,R? dz ^ c J \Pl? dz 



< 



Q(zo,R) 



< c 



Q(z ,R) 



z , R)$(u; z , R) + MR 2+2 ^ 



(3.8) 



The second counterpart of the pressure p 2 is a harmonic function and thus 
satisfies the estimate 



\P2 ~ [P 



'2\x ,Q 



2 (l.r < c,{ 



\P2 ~ [P2}x Q ,R\ 2dx 



B(x ,g) 



B(x ,R) 



\p- [p]x ,r\ dx + c J \pi\ dx 

B(x ,R) B(x ,R) 

for any < g < R. Hence, by (13. 8p . we show 
D(p 2 ;z ,g) <c[^yD(p,z ,R) + ^(u;z ,RMu;z Q ,R) + MR 2+2 ^. (3.9) 

Taking into account simple inequality 

D(p; z , g) < 2D(p 1 ; z , g) + 2D(p 2 ; z , g) < 

J \p 1 \ 2 dz + 2D(p 2 ;z ,g) 



< 



Q(z ,R) 



we deduce the estimate (13.71) from (13. 8p and (13.91) . Lemma [3^1 is proved. 

Now, we pass to the proof of Theorem 12.11 Assuming that Q(zq, R) C Q 
and r G (0, 1), we find from (I3.4f) and (13.71) two inequalities: 

z , t 2 R) < c(r 4 + z Q , tR))$(u; z q , tR) + cD(p; z , tR)+ 



+cM{tR) 2+2 ~< 



and 



[1 + c)D(p; z , tR) < (1 + c)c{t*D( P ; bz , R) + z , R)$(u; z , R) + 



+mr 2+2 ^. 

Adding the latter inequalities and introducing the new functional 

G(zo, R) = $(u; z , tR) + D(p; z , R), 

we arrive at the basic estimate 

O(z , tR) < c(r 4 + zq, tR))G(z Q: R) + 

+ c^(u; z , R)$(u; z , R) + cMR 2+2j . (3.10) 
It is not difficult to check validity of the following inequality 

zq, tR) < c$(w; z q , R) 

for any r G (0, 1) and any R > 0. Then from (13.101) it follows that 

9(z 0) tR) < c(t 4 + zq, R))Q(zq, R/t) + cMR 2+2 ~< (3.11) 

under assumption that Q(z , R/t) C Q. Letting 7x = (1 +j)/2 and choosing 
t = t(j) G (0, 1) so that 



we can choose i? < rmin{dist(f2i, cftl), vT — T{\ so that 

^(u;2 0)J R) < r 4 (3.13) 

for all zq G Qi and all < R < Rq. It is here that the modulus of continuity 
of u i— > J |w| 4 <ir<it is used. So, summarizing all the above, we have 

0(zb, ri2) < r 2+27l 6(2 , i2/r) + cMR 2+2 ~< (3.14) 

for all 2 ^ Qi and all < R < Rq. To reduce (13.141) to a known iterative 
procedure, we let g = R/r and $ = r 2 . As a result, we find 

e(z , #q) < ^ 1+7i e(z , e) + cM?? 1+ V +27 (3.15) 

for all z £ Qi an d all < p < Rq/t. The inequality (I3.15P can be easily 
iterated, see pQ, 



0(zb, ^ fc i? /r) < ^ fe(1+7) fe(z , #o/ r) + Cl Mi? 2 



h2 7 



for any k G N. Here and in what follows, all positive constants depending on 
7 only are denoted by c\. The latter inequality implies 

$(«;;8o,0*Ro) < $ k(1+ ^(e(z ,R /T)+ Cl MR 2 +2r 

for any k € N and, hence, 

$( M ;zo,^)<f? 1+7 ^ (3.16) 

for any < g < Rq, where 

\e(z ,R /r)+MRl +2 ' y 



1 \l+7 

H = Cl — 



i? 

Obviously, if is a function of i?o, IMUq, Iblb.Qj Af, dist(f2i, <9Q), T — T\, 
and 7. 

Now, our next step is to figure out how does ty(u; z ,R) depend on R. 
By (13.161) . we have 

\(uU e/2 - (u) Z0 J < -&'\u- zo, g) < cH^g-^ 2 

Q 



for any < g < Rq. Therefore, 



Rq , -(l- 7 )/2^ , 1 n(1- 7 )/2 TTl / 2 rRo\-^y 2 



for any ^ G N, or 



i=o 



!(«)«,,, - («k,*ol < CyHW-^p (3.17) 



for any < g < Rq. Proceeding and making use of (13.171) . we find 

^/ 2 (u;z ,g) < ^ 2 (u;z ,g)+cg\(u) Z0 J < c 1 g^' 2 H 1 ' 2 + 

+cg\(u) ZOiRo \. 



The latter implies 

*(u;z ,q) < Q^Hy (3.18) 

for any < g < R®, where Hi depends on the same arguments as H. 

Coming back to the basic estimate (13.10!) and taking into account ( 13. 12ft . 
( 13TT3D . (13161) . and (JHUHJ), 

6(^0, tR) < r 2+27 6(z , R) + c{M + HH 1 )R 2+2 ' y 
for any < R < R . After iterations of the latter inequality, we show 

O(z , r k R) < r k{ - 2+2 ^ (e{z , R ) + c(M + HH^R 2 ^ 
for any fcGff. Consequently, we obtain: 

Hu;z ,g)<g 2+2 m 2 (3.19) 

for any z$ G Q\ and for any < g < Rq with H2 depending on the same 
arguments as H and H\. Finally, the Holder continuity of u subject to (13.191) 
follows from known considerations, see, for example, [3] or jl]. Theorem 12.11 
is proved. 

Proof of Proposition 12.41 Our first remark is that 

L q (Qr) C M 2 , 7 (Q T ) (3.20) 

if 7 = 1 . The second remark is that 

u G L 4 (Q T ) (3.21) 

and the corresponding norm is bounded by ]| a j | 2,m 2 + ||-^1|2,q t - 

In order to handle the modulus of continuity of the function to 1— > J \u\ 4 dz, 

we are going to show that, for any < s < T, 

sup |K-,t)|| 4 ,R2 < C(s, \\F\\ 4iQt , ||a|| 2) K2). (3.22) 

s<t<T 

Assume that (13. 22[) has been already proved. Let the number r 6 (0, 1) be 
defined by (I3.12p . Then by f!3.22j) we can find a number < R\ < \yfs such 
that 

( / \u\ 4 dz J 2 < R\ sup ||w(-, t) || 4jR 2 < 

^ J ' s/2<t<T 
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< R x C{s/2, \\F\\ 4)Qt , ||o|| 2iR2 ) < r 4 (3.23) 
for any zo = (xo,to) such that t > s. Obviously, R\ depends only on s, 
\\F\\g } Q T , \\F\\2 t Q T , and ]| a j| 2,m 2 as , by interpolation, ||F||4,q t is estimated by 
||F||g Q T and ||-F||2,q t - Then we repeat the proof of Theorem 12.11 replacing Q 
with Q(z ,Ri) and Qi with Q(z , Ri/2) and establish u G L oc (Q(z , Ri/2)) 
with a uniform estimate with respect to z having t Q > s. Thus, the solution 
is Holder continuous. 

So, let us prove (13.221) . To this end, we test the Navier-Stokes by \u\ 2 u, 
as a result we have 

~d t J \u\ 4 dx + J (\u\ 2 \Vu\ 2 dx + 2\u\ 2 \V\u\f^dx = 

R 2 R 2 

= J (2\u\pu ■ V\u\ + F : (|m| 2 V-u + 2\u\u (g> V\u\)^dx. 

R 2 

After application of the Cauchy inequality with a suitable weight, we find 
the following estimate 

d t \\u\\l R2 <c / (p 2 \u\ 2 + \F\ 2 \u\ 2 ^jdx < 



kill] 



< c(||p||4 iR 2 + ll-FH^j 

It remains to make use the pressure equation which, in the case of the Cauchy 
problem, gives us 

d t \H;t)\\U* < c (ll«(^)ll8^ + ||F(-,t)||^ 2 )||n(-,t)||^ 2 . (3.24) 

To evaluate the right hand side of (13.241) . we are going to use a particular 
case of the generalized Ladyzhenskaya inequality 

||w(-,*)|li,R2 < c||u(-,*)||4 )R a||Vu||2^a. 

The proof of the generalized Ladyzhenskaya inequality is given in the Ap- 
pendix. But then (I3.24p can be reduced to the form 

1 4 



d t \H;t)\\ix> < c{\H-MU*\\vu\\i M2 + \\F(;t)\\X R2 + \\u(;t)\\b 

Multiplying the last inequality by a suitable cut-off function in t, keeping in 
mind that our solution has the finite energy bounded by ||ci||| R3 + H-F'HIqj, 
from above, and using Gronwall's lemma, we prove ( 13.221) . Proposition 12.41 
is proved. 
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4 Appendix: Generalized Ladyzhenskaya In- 
equality 

The inequality 

\\u\\l R 2 < \/2|M| 2iR2 ||V«|| 2iR2 , Vm G C °°(M 2 ), (4.1) 

was used by Ladyzhenskaya in [5J to prove unique global solvability initial 
boundary value problem for the Navier-Stokes equations in bounded domains 
of IR 2 . The generalized version of the Ladyzhenskaya inequality is as follows: 

Nll,*> < -^ll«llr I R?l|Vtt|| 2)Ra> Vn e C °°(R 2 ) (4.2) 

for r > 2. The proof of ( 14.21) is essentially the same as the proof of ( 14. If) . 
The main ingredient of it the following identity 

\u\ r (xi,X2) = r J \u\ r ~ 1 (t,X2)(\u\),i(t,x 2 )dtdx 2 - 

— oo 

Using the same identity with respect to x 2 , we find 

\u\ 2r dxidx2 < r 2 / |-u| r— |<ixi<ia;2 / \u\ r ~ \u^\dxid,X2 < 



< — J \u\ 2 ^ r r> dx\dx 2 J | Vu\ 2 dxidx 2 - 

R 2 K 2 

By interpolation, 

/ \u\ 2 ( r ~ 1 ^dxidx2 — J \ u \ 2r dxidx2 j ^ |«| r dxidx2 

K 2 R 2 R 2 

Now we deduce (14.21) from the latter inequalities. 
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